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The statistical mechanics of a mixed gas of adjoint and fundamental representation 
charges interacting via 1+1-dimensional U{N) gauge fields is investigated. In the limit 
of large N we show that there is a first order deconfining phase transition for low densities 
of fundamental charges. As the density of fundamental charges becomes comparable to 
the adjoint charge density the phase transition becomes a third order one. 


1 Introduction 

The classical Coulomb gas is an important model in statistical mechanics. It is 
exactly solvable in one dimension. In two dimensions it exhibits the Berezinsky- 
Kosterlitz-Thouless phase transition which is the prototype of all phase transitions 
in two-dimensional systems which have t/(l) symmetry. In this paper, we shall 
discuss a generalization of the classical Coulomb gas to a system of quarks which 
interact with each other through non-Abelian electric fields. This model is known 
to be exactly solvable in some special cases, Car SU{2) gauge group and fundamen¬ 
tal representation quarks in one dimension EJ and for SU{N) gaugej-group in the 
large-A limit with adjoint representation quarks in one dimension ltEI. It can also 
be formulated on thej-Lattice and solved with adjoint quarks in the large N limit 
in higher dimensions au where it has a substantially more complicated structure, 
although, even there, a solution of some special cases of the model are relevai^ to 
the deconfinement transition of three and four-dimensional Yang-Mills theory □. 

In the present paper, we shall concentrate on solving a more general version 
of the one-dimensional model than has previously been considered and elaborating 
on the properties of the solution. The one-dimensional case has the advantage of 
being directly related to a continuum field theory, the heavy quark limit of l-fl- 
dimensional quantum chromodynamics (QCD). Part of the motivation for this work 
is to study the possibility of a confinement-deconfinement phase transition at high 
temperature or density in a field theory which has some of the features of QCD, 
i.e. with similar gauge symmetries and where interactions are mediated by non- 
Abelian gauge fields. QCD exhibits confinement at low temperature and density 
with elementary excitations being color neutral particles - mesons and baryons. 
On the other hand, at high temperature or high density it is very plausible that 
the dynamical degrees of freedom would be quarks and gluons - which would form 
a quark-gluon plasma, rather than the mesons and baryons of low temperature 
nuclear physics. At some intermediate temperature or density there should be 

Talk presented by G.W.S. at the Conference on Recent Developments in Non-Perturbative 
Quantum Field Theory (AIJIC 97), May 26-30 1997 at Seoul, South Korea. 


1 



a crossover between these two regimes. There are few explicitly solvable models 
where this behavior can be studied directly. Previous to the model of H’B, known 
explicit examples of phase transitions in Yang-Mills theory were not associated with 
confinement, but were either lattice artifacts Q unrelated to the continuum gauge 
theory or werejBssociated witl^topological degrees of freedom in Yang-Mills theory 
on the sphere 9 and cylinder Qd and demark a finite range of coupling constants 
within which the gauge theory resembles a string theory. 

In finite temperature Yang-Mills theory (or QCD with only adjoint quarks), 
confinement is thought to be governed by the realization of a global symmetry which 
is related to the center of the gauge group and implemented by certain tonalpgkaLly 
non-trivial gauge transformations that appear only at hnite temperature EJ0’E3E3. 
The Polyakov loop operator is an order parameter for spontaneous breaking of this 
center symmetry and yields a mathematical way of distinguishing the confining and 
deconfined phases. When fundamental representation quarks are present, the center 
symmetry is broken explicitly and the Polyakov loop operator is no longer a good 
order parameter for confinement. Whether, in this case, a mathematical distinction 
of confined and deconfined phases exists, and indeed whether there is a distinct 
phase transition at all, is an open question. 

Here, we will consider a toy model which resembles two-dimensional QCD with 
heavy adjoint and fundamental representation quarks. It could also be thought of 
as the heavy quark limit of dimensionally reduced higher dimensional QCD where 
the adjoint particles are the gluons of the compactified dimension which get a mass 
(similar to a Debye mass) from the dimensional reduction, and the fundamental 
representation particles are the quarks. We will solve this model explmitly in the 
large-limit. The model with only adjoint quarks was solved in refs. d’H and it was 
found that the explicit solution has a first order phase transition between conhn- 
ing and deconhning phases. These phases could be distinguished by the expectation 
value of the Polyakov loop operator, which provided an order parameter for confine¬ 
ment in that case. In this paper, we shall add fundamental representation quarks. 
Then, as in QCD, the center symmetry is explicitly broken and the Polyakov loop is 
always non-vanishing. We nevertheless hnd that the hrst order phase transition per¬ 
sists when the density of fundamental quarks is sufficiently small. When the density 
of fundamental quarks is increased until it is comparable to the density of adjoint 
quarks, the phase transition becomes a second order one. When the fundamental 
quark density is increased further, the phase transition is third order. 

Another motivation of the.^resent paper, ^ well as Ref. i’i is to study a sugges¬ 
tion by Dailey and Klebanov and KutasovEfI that 1-|-1-dimensional adjoint QCD 
would be the simplest gauge theory model which exhibits some of the stringy fea¬ 
tures of a confining gauge theory. It is a long-standing conjecture that the confining 
phase of a gauge theory can be described by a string theory ll3. There are only two 
case^wj^ejthis relationship is well understood, two-diuimsional Yang-Mills the¬ 
ory EEIEZUm’Ej’O and compact quantum electrodynamics lIIlj. At low temperatures 
1-|-1-dimensional adjoint QCD is conhning in the conventional sense that quarks 
only appear in the spectrum in color neutral bound states. This is a result of the 
fact that, in one dimension, the gluon field has no propagating degrees of freedom 
and therefore it cannot form a color singlet bound state with an adjoint quark. As 
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a result, the quantum states are color singlet bound states of two or more adjoint 
quarks. The spectrum contains an infinite number of families of mnhi-niiark bound 
states which resemble asymptotically linear Regge trajectories Eil’L3’L3c3|^d, for 
large energies, the density of slate^ncreases exponentially with energy This 
implies a Hagedorn transition at high temperaturCj—Kutasov Ei supported 

this view by using an argument originally due to Polchinski that a deconfinement 
transition occurs when certain winding modes become tachyonic at high tempera¬ 
ture. This behaviour was a feature of tha explicit first order deconfinement tran¬ 
sition found in the model considered in Era. That model, which coincides with the 
heavy quark limit of adjoint QCD, is effectively a statistical mechanical model for 
strings of electric flux, with quarks attached to their ends. 

The large-expansion of two-dimensional adjoint QCD has the same complex¬ 
ity as the large-expansion of a higher dimensional Yang^ills theory and the lead¬ 
ing order, infinite-A^ limit cannot be found analytically Ell. In fact, the dimensional 
reduction of three or four-dimensional Yang Mills theory produces two-dimensional 
QCD with massless adjoint scalar quarks, so the combinatorics of planar diagrams 
is very similar. 

1.1 Overview 

This paper is organized as follows: In Section 2.1 we identify the gauged principal 
chiral model which corresponds to the gas with sources in various representations. 
In Section 2.2 the quantum mechanical formulation of this unitary matrix model is 
analyzed. This is followed by a section (2.3) where we rewrite the model in terms 
of collective variables (eigenvalue density), which are convenient for analyzing the 
large-Y limit. 

Parametrized solutions to the collective field equations are given in Section 3.1, 
and the parametrized free energy and its derivatives are obtained in 3.2. We also 
show that they give rise to a first order differential equation the free energy has to 
obey. In Section 3.3 we establish the existence of a third order line in the phase 
diagram, and compute the point where it terminates. Using numerical techniques 
we show in 3.4 that the critical line continues from that point as a first order line. 
In Section 3.5 we discuss the phase diagram we obtained. The paper ends with a 
summary (Section 4). 


2 Formalism 


2.1 Effective action 

The partition function of 1-|-1-dimensional Yang-Mills theory at temperature T 
and coupled to a number K of non-dynamical quarks at positions Xi..., Xk in 
representations ,..., Rk of the gauge group is obtained by taking the thermal 
average of an ensemble of Polyakov loop operators 


^[T, Xi,..., X}i[ , ,..., R ]^] 


dA,, 


=-S[A] 


J^TrUeMo 


(LtAq {T,Xi) 


( 2 . 1 ) 
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where the Euclidean action is 


S[A] = j dr J dx-^TT{F^i,{T,x)f , ( 2 . 2 ) 

the gauge fields are Hermitean NxN matrix valued vector fields which have periodic 
boundary conditions in imaginary time 

x) = A^{t + 1/T, x) , 


and the field strength is 


Ffiu — dyA^ i [j4^, Ay^ . 

The gauge field can be expanded in basis elements of the Lie algebra A^' 
with the generators in the representation Ri. For concreteness, we consider 
U{N) gauge theory and denote the generators in the fundamental representation as 
T“ with a = 1,..., iV^. They obey 





(2.3) 

normalized so that 

Tr 

= 2^ ■ 

(2.4) 

and with the sum rule 

A72 

II 

to 1 



\ rjiCL rpQj 

/ 

a—1 

(2.5) 


We remark that group elements in the adjoint representation are related to the 
fundamental representation matrices g by 


(^Ad)ah ^ 2 Tr(gtT“gr'’) . 


( 2 . 6 ) 


The expression ( |2.l| ) can be obtained by canonical quantization of 1+1-dimensional 
Yang-Mills theory with Minkowski space action coupled to some non-dynamical 


sources 


S=- 


dtdx —TrF^yF^’' 


+ source terms . 


(2.7) 



In the following, we will review an argument for representing the partition func 
of Yang-Mills theory as a gauged principal chiral model which was first given ir 
and which was generalized to the case of Yang-Mills theory with sources in I 
is usual in canonical quantization of a gauge theory, the canonical conjugate E(x) of 
the spatial component of the gauge field (which we denote by ^(x)), is proportional 
to the electric field, 

1 

'+1 


E = -y+n 


and obeys the commutation relation 


[A%x),E\y)\ = z(5“'’5(x - 2 /) 
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( 2 . 8 ) 



The Hamiltonian is 




/ dx-Y,{E%x)f . 


(2.9) 


a=l 


This Hamiltonian must be supplemented by the Gauss’ law constraint equation 
which is the equation of motion for Aq following from (2.7) and which contains the 
color charge densities of the sources 


K 


g-(x) = ( ^E^x) - r^^A\x)E-ix) + ^ TlSix - x.) ) - 0 . 


( 2 . 10 ) 


Here, the particles with color charges are located at positions Xi,... ,xk- are 
generators in the representation Ri operating on the color degrees of freedom of the 
i’th particle. 

There are two options for imposing this constraint. The first is to impose 
another gauge fixing condition such as 


H ~ 0, 


and to use the constraints to eliminate both E and A. The resulting Hamiltonian 
is 

g2 JY 

H=J2 ® - ( 2 - 11 ) 

i<j,a 

which was considered in Ref. 0. It is the energy of an infinite range spin model 
where the spins take values in the Lie algebra of U{N). 

The other option, which makes the closest contact with string dynamics, is to 
impose the constraint (^.IC ) as a physical state condition. 


g‘'(x) ^'phys = 0 . 


To do this, it is most illuminating to work in the functional Schrodinger picture, 
where the states are functionals of the gauge field, iplA] and the electric field is the 
functional derivative operator 

The time-independent functional Schrodinger equation is 


/ 


dx 



5^ 


(<5H“(x))2 






Gauss’ law implies that the physical states, i.e. those which obey the gauge con¬ 
straint ( 2.10D , transform as 

[A^;xi,...,xk] = g^^\^{xi)...g^^^^{xK)'i’‘’^-''^[A;xi,...,XK] , 
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where 

= gAg^ - ig'^g^ , 

is the gauge transform of A. 

For a fixed number of particles, the quantum mechanical problem is exactly 
solvable. For example, the wavefunction of a fundamental representation quark- 
antiquark pair is 

= Ipe'-f-TV’ , (2.12) 

where the path ordered phase operator represents a string of electric flux connecting 
the positions of the quark and anti-quark. The energy is — X 2 \- For a pair 

of adjoint quarks, the wavefunction is 

■f‘^^[A;xi,X2] = Tr fT‘^re -^-1 . (2.13) 


with energy ^-^\x 2 — xi\. These energy states are identical to what wo uld b e 


obtained by diagonalizing the ‘spin’ operators in the gauge fixed Hamiltonian (2.11). 

Note that the wavefunctions ( 2.12| ) and ( 2.13 ) are not normalizable by functional 
integration over A. This is a result of the fact that the gauge freedom has not been 
entirely fixed, so that the normalization integral still contains the infinite factor of 
the volume of the group of static gauge transformations. 

In general, for a fixed distribution of quarks, a state-vector is constructed by 
connecting them with appropriate numbers of strings of electric flux so that the 
state is gauge invariant. The number of ways of doing this fixes the dimension of the 
quantum Hilbert space. If the flux strings overlap, the Hamiltonian can mix different 
configurations, so the energy eigenstates are superpositions of string configurations. 
However, this mixing is suppressed in the large-fV limit (i.e. the strings are non¬ 
interacting) and any string distribution is an eigenstate of the Hamiltonian with 
eigenvalue (e^A^/4)x(total length of all strings). 

We shall study the thermodynamics of this system by constructing the partition 
function. We work with the grand canonical ensemble and assume that the quarks 
obey Maxwell-Boltzman statistics. The partition function of a fixed number of 
quarks is constructed by taking the trace of the Gibbs density over physical 

states. This can be i mpl emented by considering set of all states in the representation 
of the commutator (pT^), spanned by, for example, the eigenstates of A°'{x) and an 
appropriate basis for the quarks 


1^) ® Boi <8) Ca 


Projection onto physical, gauge invariant states is done by gauge transforming the 
state at one side of the trace and then integrating over all gauge-|transformations 
(and then dividing by the infinite volume of the gauge group) S The resulting 
partition function is 

Z[T/e‘^;xi,...,XK] = J [dA][dg] {A\ \A^) Tr (xi).. .Tr (xk) , 

(2.14) 
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where [dg\ is the Haar measure on the space of mappings from the line to the 
group manifold and [dA] is a me asure on the convex E ucli dean space of gauge field 
configurations. The expression ( 2.14 ) is identical to (|2.l|) with the Polyakov loop 
operator is the trace of the group element g{x) in the appropriate representation. 

In going over to the grand canonical ensemble the first step is to integrate over 
all particle positions. We then multiply by the fugacities for each type of charge: 
a factor of A/j for each quark in representation R. To impose Maxwell-Boltzmann 
statistics, we divide by the factorial of the number of quarks in each representation. 
We then sum over all numbers of quarks in each representation. This exponentiates 
the fugacities, resulting in the grand partition function 


Z[T/e^\R] = 


J [dA][dg] 




(2.15) 


where the effective action is 


exp = {A\ e |A®) exp dx XrTt g^{x)^ , (2.16) 

and the summation in the exponent is over all the irreducible representations of 
U{N) we want to consider. The Hamiltonian is the Laplacian on the space of gauge 
fields. The heat kernel obeys the equation 

with the boundary condition 

^^hm ^ (H| |H®) = (5(H - H®) . 

These equations are easily solved by a Gaussian - divided by a T-dependent con¬ 
stant: 

{A\ e-^/^ \A9) ~ exp (^- J dx ^ Tr {A - H®)^^ . 

We see that the effective theory is the gauged principal chiral model with a 
potential energy term for the group-valued degrees of freedom. 






Vg{x)+i[A{x),g{x)] {XrTt g^{x)) 


R 


Note that we have introduced the coupling constant 

N 


7 = 


2T 


(2.17) 


(2.18) 


When we analyze the limit N —> cx), we will tune such that 7 is constant. 
Moreov er, w e assume that the fugacities Xr are scaled such that all terms in the 
action ( 2.17 ) are of order . 
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The potential energy term in the effective action, 


V{g) ^ -^A«Tr(g«(x)) , (2.19) 

R 

is the expansion of a local class function of the group element g{x) (one which 
obeys V{g) = V{hgh~^) for h G U{N)) in group characters with coefficients A^;. 
The characters 

Xnig) = TT{g^{x)) , 

form a complete set of orthonormal class functions of the group variable, with inner 
product 

J [dg]XR{9)XR'{9) = Sr.r' . 

Here [dg] is not a functional integral measure, but is the Haar measure for integration 
on U{N). From the potential, we can find a fugacity by 

= - J[dg]x*R{9)y{9) ■ 


By tuning the fugacities appropriately, we could obtain any local invariant potenli^ 
The effective action (2.17) with all A_r = 0 was discussed bv£rignani et.al. Ej’lZI 

00 by ^rembo The model with 


and was solved explicitly in the limit N 
XacL ^ 0 (with adjoint quarks) was solved in Refs. □ and cl. The effective action 
( 2.17 ) is gauge invariant, 




It is also covariant under the global transformation 

Ses[Az9,XR] = Ses[A,9,z-^dR)Xj,]. ( 2 . 20 ) 

where 2 ; is a constant element from the center of the gauge group, which for U(N) is 
^ U{1) and would be the discrete group Zn for gauge group SU{N). Here, Ci{R) 
is the linear Casimir invariant of the representation R, which is the number of boxes 
in the Young tableau corresponding to R. When the gauge group is SU{N) and 
the only non-zero fugacities are for the zero ‘N-ality’ representations, i.e. those 
for which Ci{R) = 0 mod N, there is a global Z^- symmetry. For gauge group 
t/(lV), this occurs only when all representations with non-zero fugacities have equal 
numbers of quarks and anti-quarks. The fugacities of other non-symmetric charges, 
can be thought of as an external field which breaks the center symmetry of the 
system explicitly. This situation is akin to the effect of an external magnetic field 
on a spin system. 


2.2 Matrix quantum mechanics 

If we re-interpret x as Euclidean time, the partition function that we have derived 
has the form of a Euclidean space representation of the partition function for matrix 
quantum mechanics, where the free energy is identical to the ground state energy 




of the matrix quantum mechanics. We can study the latter model by mapping the 
problem to real time r by setting x = ir and A —> —iA. The action in real time is 
then 


Sqm 



^T,\g + ^[A,gf 



We remark that this action must not be confused with the action (2.2). Sqm is 
the action for a 0+1-dimensional problem (quantum mechanics), while (2.2) is the 
action for Yang Mills theory in 1+1 dimensions. This remark also holds for the 
Hamiltonian below. In order to avoid confusion, we label the quantum mechanical 
quantities with the subscript QM. 

The canonical momentum conjugate to the group valued position variable g is 
the Lie algebra element 


n = — + a'^Ag - A) , 

and the Hamiltonian is 

Hqm = ^Trtf - ^TrUig^Ag - A) + V{g) . (2.21) 


The gauge field A plays the role of a Lagrange multiplier which enforces the con¬ 
straint 

gAg'^ - H ~ 0 , 

and the Hamiltonian reduces to 

Hqm = :^Trn2 + V{g) . 


We can expand the canonical momentum as 


H = ^ n“r“, 

a 

Then, the components satisfy the Lie algebra 


n“,n^] 

= , 

(2.22) 

[n“,5] 

= aT\ 

(2.23) 


= -T‘^a^ . 

(2.24) 


It follows that in the Schrodinger picture the components of the canonical momen¬ 
tum are represented as 

n- = . 

Denoted in components the constraint reads 

G“ = ^TrT“gT^g^ - H^ - 0 . 


9 









The constraint has no operator ordering ambiguity. It generates the adjoint action 
of the symmetry group 

[G“,5] = • 

The constraint can be realized as a physical state condition 

V’phys = 0 . 

In the representation where states are functions of g, this implies that the physical 
states are class functions 

'*Aphys(5') “ '4’physihgh ) , 

where h € U(N). This means that the physical states are functions of the eigenval¬ 
ues of g. In a basis where g is diagonal, 

/e“i 0 0 0 ... 0 \ 

0 0 0 ... 0 

g = 0 0 0 ... 0 , (2.25) 

0 0 ... 0 

Vo 0 ... 

the wavefunctions are functions of Oj, 

'^phys — V^phys (o^l i ■ ■ • j “t“ 27r, Ojv). (2.26) 

Denoting the gauge group Laplacian in components 

A = ^(n“)2, (2.27) 

a—1 

the Hamiltonian reads 

Hqm = -h V{g) . 

Since the potential V(g) is also a class function and depends only on the eigenvalues, 
when operating on the physical states, the Hamiltionian can be expressed in terms 
of eigenvalues and derivatives by eigenvalues 

= 3^ jV) (? ^ • 

where 

T/ \ TT n • \ ^ e/(Q!) 

'Ho:) = jy[2sm-(Q;i - aj) = (2j)3V(3V-i)/2 ^(Af-i )/2 > 

and 

j{z) =W_Hi - ^i) ) Zi = e*“A 

i<j 
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The physical fStales must by symmetric functions of ai. (There is a residual 
gauge invariance E3o under the Weyl group which permutes the eigenvalues and 
the physical state condition requires that the physical states be symmetric under 
these permutations.) The normalization integral for the wavefunction is 

J = 1 ■ 

Since the integrand depends only on the eigenvalues of g, It is convenient to write 
the Haar measure as an integral over eigenvalues of g with the Jacobian which is 
the Vandermonde determinant, 

f {Y[dai)\J{a)ftp''{a)ilj{a) = 1 . 

i 

The Hamiltonian and inner product have a particularly simple form when we rede¬ 
fine the wavefunction as 

... ,aAr) = J(a)V^(ai,.. .. un ) ■ 

Since J is antisymmetric, '0 is a completely antisymmetric function of the eigenval¬ 
ues, which we can think of as the coordinates of fermions. The Hamiltonian is that 
of an interacting Fermi gas 

|i^ (e + ^(«)| ^(«) = • 

This cop:espondence of a c=l matrix model with a Fermi gas was first pointed out 
in Ref. E. 


2.3 Large N: Collective variables 

In this section we shall examine the collective field forui iilatir i Ti p f the large-TV limit 
of the theory that we discussed in the last subsection □O’EE. The Hamiltonian 
obtained in the last subsection reads 


H, 


QM 


_7 

47V 




+ V{g) 


(2.28) 


a=l 


with (compare (^.I9|)) 


V{g) = -^AflTr( 5 «(r)) . 


R 


It was shown (compare ( 2.25 ),( 2!^ )) that the wavefunction depends only on the 
eigenvalues of g and thus the density of eigenvalues 


1 

= — ^J(6»-a,(r)) , 


11 












completely characterizes the properties of the system. Interpretation of the physics 
of the system at large N is more convenient when one considers the Fourier trans¬ 
form of the eigenvalue distribution 

n/0 

where we have defined 

Cn(r) = ^Trg”(r) , c-nir) = Cnir) . 


We now turn our attention to developing the collective field theory formulation 
of the (thermo-) dynamical problem given by the Hamiltonian ( ^.28 ). Since the 
wavefunction depends only on the eigenvalues of g, we would like a Hamiltonian 
equivalent to ( |2.28| ) but written in terms of the eigenvalue density p and a conjugate 
momentum H. At large N we will find this Hamiltonian and write equations of 
motion for p and H. So far we have not imposed any restriction on the potential 
V{g)^ but from now on we assume, that it can be expressed as a functional of the 
eigenvalue density p{9). In particular the potential we are going to analyze below 
will have this property. 

The canonical momentum operates on the wavefunction as 


n“V’[p]= / d0[n“,p(0)] 


6p{6) 






and the Laplacian (2.27) is 
1 


AV’[p] = 


K 


K 


I 47riV J 

1 


\K\ Tr/Trc/^"-^ - NTig^ 


K 


\L=0 


d0d0' 


KL 


Sp{0)5pi0') 


5p{0) 

V'[p], 


which can be written as 


= -^ J d0p{0)^ 


d 5 
m 5p{0) 
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where 


m = ^'P I de'p{9') cot (^) • 

V indicates principal value integral. 

The transformation of the wavefunction 




bb]exp 



dOdO' In sin 


\e-e'\ 


p{9)p{0') 


transforms the derivative in the Schrodinger equation so that it has the form 

+ *1 = eH,]. (2.30) 

The second term on the left-hand-side of this equation has a simple form. In 0 
it is shown that it gives rise to a term which is cubic in the density. Thus, the 
Schrodinger equation has the form 

(2.31) 

up to an overall constant. 

The large-fV limit is dominated by the eikonal approximation. In this approxi¬ 
mation, we make the ansatz for the wavefunction 

tp[p] = exp (ifV^S'b]) . 


The eikonal, S then obeys the equation 


Hqm[p,^] 

7V2 


7 

8 


d9 { p{9) 


f d 5 S 
\d9 5p{9) 





E 


(2.32) 


Here, we have ignored a term which is of subleading order in N'^. We have also 
assumed that V\p\ will be of order (compare Section 2.1) and that the natural 
magnitude of the energy eigenvalue is of order . 

To solve this equation for the ground state, we must find its minimum by varying 
p and the canonical momentum 


n = 5S/5p, 


subject to the condition that p is normalized. This leads to the equations of collec¬ 
tive field theory 


dr 


p{t, 9) 


SHqm/N^ 
(5n(r, 9) 


^v(r,e) 


^ 6Hqm/N^ 
d9 Sp{t, 9) 
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where 


e) = -^ n(r, 6). 

Taking the derivative of the second equation with respect to 9 eliminates a Lagrange 
multiplier which must be introduced on order to enforce the normalization condition 
for p. Using (p.32|) one finds 


dr 


dr 


7 d 




7 d 

I'm 

_ 

Ir^mTp 


1 


{pv) 

V[P\ 


0 , 

0 . 


(2.33) 


It is interesting to note that these are nothing but Euler’s equations for a fluid with 
equation of state P = p^ jZ on a cylinder with coordinates (0, r). The inclusion of 
a potential V(9, r) corresponding to non-Abelian charges is equivalent to subjecting 
the fluid to an external force which is derived from V(9,t). 

We shall use these equations in the next section where we analyze the large- 
limit of a mixed gas of adjoint and fundamental charges. 


3 Free energy and critical behaviour 

3.J Static solutions to the collective field equations 


In this section we will find static solutions to the collective field equations (2.33). 
The most simple potentials involve only the lowest representations, the fundamental, 
its conjugate and the adjoint. We shall consider a slight generalization of these and 
use powers of the lowest representations to include multiple windings of the Polyakov 
loop operator. Our potential reads 

OO 

V{g) = - ^ (K„iVTr(g") + k„iVTr((gt)") + A„|Trg"|2) , (3.1) 


where we made use of (compare (|2.5|), (2.6)) 

Tr {g^<^{x)r = |Tr g^x)f , 

to relate the trace in the adjoint representation to the trace in the fundamental 
representation. The couplings for the fundamental representation charges (and 
their conjugates) were chosen to scale ~ N, to make the potential of order N'^. It 
should be remarked that parts of this section can easily be extended to more general 
potentials. 

The potential o indeed can be expressed as a functional of the eigenvalue 
density (2.29|). The collective field Hamiltonian (^.32|) then reads 




QM 


7V2 


= Z d9 


Pi0) + 


^ ( ^r, 


d9 p(9) 


AnO 


+ Kn d9 p{9) e“^ + Kn d9 p{9) e 


— inO 


(3.2) 


96 
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In order to maintain correspondence with the original version of the Hamiltonian 


(2.28) we subtract the constant 7 / 96 . It sets the energy scale such that the free 
energy vanishes in the confined phase of the model with only adjoint charges (see 
below). 


The corresponding collective field equations (2.33) read 


dp , 1 d , X n 


(3.3) 


dv 7 d ^ r/> 


+ Kn) e™^ + {XnCn + Kn) 6 — 0 . 


(3.4) 
Note 
iv in these 


where c„ are the x-dependent Fourier coefficients of p as introduced in 
that we also performed the change of variables, r —> —ix and v —■ 
equations in order to invert the Wick rotation performed at the beginning of Section 
2.2 prior to canonical quantization. 


We will only consider real, static solutions of the non-linear equations (3.4), 
that is, where p{0,x) = Pq(9) and the velocity v vanishes identically. Consequently 


+ Y.i^nC-n + Kn)e'^^^ + J^i^nCn + K„)e-“® where p is real 
0 otherwise 

(3.5) 

The constant of integration E has physical interpretation as the Fermi energy of a 
collection of N fermions in the potential V[p\ and is fixed by the normalization 
condition 

1 = y po{e) . (3.6) 

Here it is more convenient to express the c„ in terms of po (compare ( 2.29| )) 

Cn = j d9 po(0)e' 


And 


(3.7) 


The real support of the function pq{Q') is the positive support of A = i?-|-^(A„c_„-|- 
+ '^{XnCn + The zeros of A define the edges of the eigenvalue 

distribution and when these zeros condense one has critical behaviour in the ob¬ 
servables of the model as in general Hermitean and unitary matrix models. 


3.2 A differential equation for the free energy 

In this subsection we compute all first order derivatives of the free energy and show 
that they obey a differential equation of the Clairaut type. 

Inserting the static solution (|3.5D in (|3.2|) we obtain for the free energy 


1 11 

s / = jE - ji: 

n—1 


2 ( 


K„c_„ ) - ■ (3.8) 
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Note that / is the leading coefficient {0{N‘^)) of the energy for the matrix quantum 
mechanics problem, but for the quark gas problem plays the role of the leading 
coefficient of the energy density. 

Deriving the expression with respect to Aj for some fixed J and using 

derivatives of Equations ( |3.6| ) and (3.7) with respect to the same parameter one 
obtains 

rif 

(3.9) 


df I ,2 

517 = -'"•'I 


We use the notation d/dXj to indicate that also the Cn and E which implicitly 
depend on Xj are derived with respect to this coupling. Similarly one can show 


A. 

dnj 


= -cj 


and 


dj 37 


00 


2 An|Cn| KfiCn KnC—r, 


1 

% 


(3.10) 


(3.11) 


It is interesting to notice, that combining Equations (^) - ( 3.1l| ) gives rise to a 
first order differential equation of the Clairaut type 


7 


d'y 


00 

n 


, df df _ df 

AnTJT-1" -1“ 

Cl\'ri diKri duK/ri - 


= /■ 


(3.12) 


This differential equation has general solutions of the form 

7 7 


/(7,A„,«„) = 7F(^,^), 


where E is some arbitrary smooth function. This result shows that the parameter 
7 is not driving the physical properties of the model, but rather sets the energy 


scale. The differential equation (3.12) gives no further restrictions on the function 
F and another analysis will be adopted in the next section. However, when find ing 
the physical interpretation of the phase diagram the differential equation ( 3.12| ) is 
a valuable tool. 


3.3 Regime of the third order phase transition 

We now restrict ourselves to the case of only one pair of non-vanishing couplings 
Xj, Kj 7^ 0. Furthermore it is sufficient to consider kj real, since an eventual phase 
of Kj can always be removed by using the covariance ( 2.20 ) of the action and Haar 
measure [dg] under tran sformations by a constant element of U{1). 

In the form of (3.5)it is evident we need to solv e sim ultaneously for the nor¬ 
malization condition (3.6) and the Fourier moment (3.7) in order to have a self- 
consistent solution of the saddle-point equations. We begin by introducing an aux¬ 
iliary complex parameter 


rj e 


iJ0j — 


= Xj C-J + Kj 


(3.13) 
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and rescaling the Fermi energy as 


E = 2^irj . (3.14) 

With this notation the normalization and moment equations are respectively 


1 = 2 


Hp) 


cje 


iJjS 


= 2 


M{fi) 


(3.15) 


where we have defined the integrals 

/(m) = 

M(/r) = 


\//i + cos{J9) H(/i + cos{J9)) dO 


f cos J9^j /r + cos{J9) H(^ + cos(J6*)) d9 . (3.16) 

J —TT 


Here H(..) denotes the step function. A simple transformation of the integration 
variable shows that J(/i) and M(/i) are independent of J. Thus J enters only as 
the subscript of the parameters. For notational convenience we abbreviate 


Xj = X 


Kje = K 


(3.17) 


We remark that /(/i), and thus are real. Eliminating the moment cj 

from ( 3.15|) by using the definition (^.13 ) we obtain 


K 

1 




T - 
4 7 


(3.18) 


This family of lines in the k, A-plane parametrized by ^ represent a necessary condi¬ 
tion which a solution of the normalization and moment equations (3.15) must obey. 
From the last equation it is obvious that also the product k = is real 

(thus is just a sign). It occurs as a natural parameter when rewriting the free 
energy in terms of /(/r) and M(^) (use (3.8)) 


/ = 


12/(/r)2 


2/r — 


M{fi) 


I { 9 ) 


Mjfi) 


96 


(3.19) 


where we have eliminated A using the necessary condition (3.18). Also the first 
derivative of the free energy with respect to k can be expressed conveniently in 
terms of I(p) and M(/r) 


dn 


= -2 


M(m) 

^(m) 


(3.20) 


Remember that we restricted ourselves to kj real, and thus we encounter a factor 
2 compared to 3.10), since a real kj is the same for both terms cj and c_j in the 
potential (3.8). 

With the parametric solution ( |3.8| ) at Jaand we turn our attention to establishing 
the critical behaviour in this model. In a it is shown that the first derivati ves o f 
I{fi) and M{ii) have non-analytic behaviour at fi = 1, hence the expression (|3.2C|) 
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suggests that the vicinity of fj, = 1 is a natural place to look for non-analytic 
beh^our in the f ree energy of our model. Using the explicit results for 7(1), M(l) 
(see E3) and ( ^.18 ) we obtain the necessary condition for the critical (/r = 1) values 
of A and k 

^ ^ (321) 

7 512 3 7 ■ ^ 

Having identified a line in the phase space where we expect critical behavipuj j 
will now proceed to establish the details of this critical behaviour. Following c3c3' 
we begin by expanding about /i = 1 and the line (^.21) 


/i — l-t-£,£^>0 , 7(l-t-e) — Iq SI 


M{l + e) = M^ + SM , (3.22) 


(7c = 7(1), Me = M(l)) and analyze the variation of k around while keeping A 
fixed 

k = k'^ + 5k , A = A^ . (3.23) 

Due to the nontrivial support of the integrands in ( 3.16| ) in principle one has to 
distinguish the cases /i > 1 and ^ < 1; (see o for details). In orde r to k e ep th e 
formulas simple, we explicitly analyze only the case ^ > 1 as given by ( |3.22 ), ( 3.23 ). 
The case /i < 1 can be treated along the same lines and we denote the corresponding 
results in the end. 

The expansion now consists of two steps. We first expand the necessary condi¬ 


tion (3.18) at /r = 1 to obtain the relation between the variation 5k, and e. In the 
second step we expand the right hand side of ( f3.2C| ) at /r = 1 and use the result 
of step one to express the variation of df/dn in terms of Sk. The latter result can 
then be used to analyze eventual singular behaviour of higher derivatives of the free 

energy. _ _ 

Expanding the necessary condition ( ^3.1^ ) and using ( 3.21 ) we obtain for the 
variation of k to lowest order 


Sk 

1 51 

M Me 

■ 5M 

51' 


■-TT^ 4A^' 

1 

2 {Ic? 

7 Ic 

[ Me 

Ic\ 


256 ' 3 7 _ 


M 

7c 


We 

7 2 


(3.24) 

where in the last step we made use of the relation between the variations SI and 
5M and inserted the explicit results for L = 7(1) and Me = M{1). Using the result 
for 51 jlc to lowest order we obtain (see a for further details) 


Sk = — e ln(£) a 

where we introduced the abbreviation 

1 ^ 


a = 


77r^ 4A= 


(3.25) 


(3.26) 


Inverting equation (3.24) (again taking into account only the leading order) gives 

^5k [in ^ . (3.27) 


e = — a 


This equation is the relation between the variation Sk and e which is implied by 
the necessary condition ([3.18|). In the final step we expand the derivative of the 
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free energy (3.2C) at /r = 1 and use the result (3.27) to obtain the variation of the 
derivative in terms of Sn. Expanding (3.20|) gives 


dn 


= - 2 - 


.M, 


1 + 


5M 


51 

Z 


= -5 - ' 


Using (3.27) we obtain 


dK 


“ "3 + 


-cr ^6k + a ^Sk [in (tr ^(5 k)] 


(3.28) 


(3.29) 


We remark, that the case ^ < 1 with expansion fj, = 1 — e, e>0 changes only the 
sign of the argument of the logarithm. Differentiating the last result with respect to 
6k establishes the singular behaviour of the third derivative of the free energy with 
respect to k. Thus we find a third order phase transition for ^ = 1. The critical 
line is a straight line given by (|3.21). 


It is important to notice, that at (see Equation (3.24)) 

W 37r2 
7 


1024 ’ 


(3.30) 


the leading term in the expression for Sk vanishes. Equation ( p.24 ) is reduced to 
the simpler relation 

- = -E ^ . (3.31) 

7 27 

At this point the expansion of df / dn gives 

I = _| _ l£l„(E) - r = -| + (3-32) 

Again the case /r < 1 differs only by the sign of the argument of the logarithm. 
Differentiation with respect to Sk show s, that the phase transition has turned to 
second order at that point. Using ( 3.21 ) one can compute also the k/j coordinate 
of the second order point giving k/j = 7r^/1024, A/7 = 37r^/1024. In fact the more 
global analysis in the next section will show, that the third order line terminates 
at the second order point k/j = 7r^/1024, A/7 = 37r^/1024, and continues as a first 
order line. 


3.4 Regime of the first order phase transition 

As pointed out at the end of the last section at the point k/'j = 7r^/1024, A/7 = 
37r^/1024, the third order transition along the /r = 1 line ( |3.2l|) changes to second 
order. This unusual behaviour requires further investigation which we will carry out 
in this section. To begin, a graphical analysis of the phase diagram is mo st us eful 
and in Figure ^ we plot a number of representatives of the family of lines ( ^.1^ ) for 
a range of values of pL. 

It is clear that in most of the k, A-plane points are in a one-to-one correspon¬ 
dence with values of the parameter fi. This correspondence breaks down though in 
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Figure 1: Plot of the lines (3.18) for ji ranging from 0.4 (upper right corner) to 75 (line at the 
extreme left). The region of overlapping lines corresponds to a region of first order phase transition. 


a small region near the A /7 axis between A /7 = 0.05 and A /7 = 0.0625. Due to 
the behaviour of the slope and intercept in the linear equation (3.21) lines begin to 
overlap for increasing n starting at /i ~ 1 and continuing as /r —*■ 00 . In this overlap 
region the phase diagram is folded at the vertex k/j = 7r^/1024, A /7 = 37r^/1024, 
and each point falls on three different lines of constant /i. Consequently the system 
simultaneously admits three configurations with different free energies in this region 
of the phase space. This circumstance allows for a first order phase transition to 
develop along a line where the free energies of the different phases are equal. 



K/y 


Figure 2: Plot of the boundary of the multiple phase region. The boundary (solid curve) is given 
by a caustic of lines in the one-parameter family (3.18). The dotted curve shows the numerically 

determined first order line. 


The edges of the triangular first order region in Figure ^ is given by a caustic of 
lines from the one parameter family (|3.18| ). The boundary is defined by the curve 
where the family of curves is stationary with respect to fj,. This condition can be 
used with (3.18) to give a definition of the boundary caustic shown in Figure 2. The 
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stationary condition can be solved with the parametric result 


K _ 1 

7 4 / 2 (^) ^ ^ 

As can be seen, the curve given by ( tl.33[ ) intersects the A/y axis at two points: 
0.057024 (/i = 0.95324) and 1/16 (/i = oo) and reaches a singular maximum in the 
k /7 direction for ^ = 1 at the point re/y = 7r^/1024. The end of this region of first 
order transitions agrees with the position of the second order transition point which 
was determined by the analysis of critical behaviour in the previous section. 

Once one has determined the region where different phases can co-exist the next 
issue to address is that of the position of the line of first order phase transitions 
where different phases have the same free energy. 



X/y 


Figure 3: Free energy //y as a function of A/y in the region of first order phase 
transitions. Each curve is plotted for fixed re/y which from right to left is re/y = 
0, 0.0005, 0.001, 0.0015, 0.002, 0.0025 


The first order line can be determined for given At/y by the simultaneous solution 
for the parameters /ii and /i 2 of the pair of equations 


4 /(/ii)/(^2) 




(3.34) 


and 

1 / Ml _ ^^2 \ 1 / M{yi2) \ _ K I{yi2)M{fii) - M{p,2)I{p-i) 

6 1 /(^l)/(/l2) 

(3.35) 

Unfortunately, these equations are analytically intractable. Again we turn to 
a graphical analysis to gain further insight. In Figure we plot the free energy of 
the system as a function of A/y for different values of fixed K/y. From here it is 
easy to see a number of features of the region of first order transitions. Increasing 
/I traverses these curves in a clock-wise rotation so that the free energy increases 
for small values of /i, intersecting the nearly horizontal large yi free energy. This 
intersection point is a graphical demonstration of the first order transition which 


21 















occurs here as the model jumps from weak (/i < 1) to strong coupling large). 
Each phase continues to exist after the transition point and may be reached by 
an adiabatic process until ending in cusps which mark the boundaries of the first 
order region in the A/y axis. It is interesting to note that there is an energetically 
unfeasible intermediate “medium coupling” phase which connects the weak and 
strong phases. Hence, for fixed At/y there exist three distinct conhgurations of the 
system for given A/y in the region of first order transitions. 

3.5 Phase diagram 

We end this section with a short discussion of the phase diagram analyzed in the 
previous sections. Figure 4 shows a schematic picture of the phase diagram. The 



Figure 4: Schematic picture of the phase diagram. The dashed curve marks the first order part 
of the critical line. The solid curves above and below it are the boundaries of the area with two 
possible phases. They join at a point which shows second order behaviour. For larger re/y we find 
a third order line (p = 1) marked by a solid line. The numbers label four extremal corners of the 
phase diagram. We also indicate the range of the auxiliary parameter p. 


dashed line gives the numerically determined hrst order line. The full curves above 
and below show the boundaries of the region where different phases can coexist. 
With increasing tt/y, this region narrows, and ends in a point which is second order 
as discussed above. The curve separating the two phases continues as the third 


order line (3.21). 

We also show the range of the auxiliary parameter fi. In the hot and dense 
(deconfining) phase it assumes values /t < 1, while in the cold and dilute (conhning) 
phase it is restricted to /r > 1. On the third order part of the critical line we have 
^ = 1 . 

It is interesting to discuss the extremal regions of the phase diagram: In the case 
of only adjoint charges (all = 0) the action is invariant under transformations in 
the center of the gauge group, g ^ zg (compare ( ^.2C| )). When this symmetry is 
faithfully represented we expect 


Cn(x) = lim 


N- 


— (Tr o" 


(x)) =0 Vn 7^ 0 


(3.36) 
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since 


Tr5"(a:) ^ z” Tr g"(x) 


From (^), ( |3.4| ) it is obvious that ( ^.36 ) which corresponds to constant eigenvalue 
density p = l/27r is a solution. The stability analysis performed in □ shows that this 
solution (’’confined phase”) is stable for A/y £ [0,1/16]. Since the upper boundary 
of the multiple phase region is at A/y = 1/16 we confirm this result. 

The center symmetry of the action is explicitly broken when there are funda¬ 
mental charges (some Kn ^ 0). In particular this is the case when we consider 
the model with only fundamental charges (and their conjugate). Thus we cannot 
expect to find a confining solution with vanishing c„ which would correspond to 
vanishing (energy) densities of colour electric string {ppiux), fundamental repx^en- 
tation charges (pp,p) and adjoint representation charges (pAd) (please see llj for 
details). For this case (A = 0) we have established the existence of a third order 


phase transition at «;/y = 7r^/512 (set /r = 1, A = 0 in (3.18)). 


Finally we discuss extremal corners of the phase diagram, which are labeled 
1,..,4 in Figure The four extremal cases can easily be understood by the magni¬ 
tude of the energy densities ppiux, PAd and pp p- Point 1 is in the extremal corner 
of the low density phase. All three densities are rather small. Points 2 and 3 are 
both in the high density phase, in areas which are dominated by adjoint charges 
(Point 2) and fundamental charges (Point 3). It is nice to see, how the energy 
density of the sources is dominated by the contributions of the adjoint charges and 
fundamental charges respectively. Finally Point 4 is in a region where both the 
density of the adjoint charges as well as the fundamental charges is high and of the 
same magnitude. 


4 Summary 

In this paper we analyzed the thermodynamic properties of a model of static sources 
on a line interacting through non-Abelian forces. It was shown that the partition 
function takes the form of the partition function of a gauged principal chiral model. 
Using the eigenvalue density as collective field variable the Hamiltonian for the 
eigenvalue density in the large A^-limit was computed. We gave a static solution 
of the corresponding Hamilton equations. For the special case of only two types of 
charges, the static solution was parametrized using the parameter p proportional 
to the Fermi energy. In particular the case of two types of charges transforming 
under the adjoint, and charges transforming under the fundamental representation 
of the gauge group was considered. Expanding the parametrized solutions at /i = 1 
we established the existence of a straight line in the phase diagram where the free 
energy exhibits a third order phase transition. We proved that the third order 
behaviour terminates at a second order point. The critical line then continues as a 
first order line, which was determined numerically. 

The whole phase diagram was interpreted by analyzing the contributions of 
charges and flux to the energy density. We found that for /r > 1 the system is 
characterized by low energy densities, while for p < 1 the densities are high. 

This 2-dimensional model could be generalized in several directions. It would be 
interesting to analyze non-static solutions of the Hamilton equations and different 
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boundary conditions which might be used to include a 0-term. Loop expansion of 
the fermion determinant of QCD 2 with large quark masses could be used to relate 
the fugacities of the non-Abelian gas analyzed in this article to the mass parameters 
of QCD 2 . 
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